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Gauge-invariant systems in unconstrained configuration and phase spaces, equivalent to second- 
class constraints systems upon a gauge-fixing, are discussed. A mathematical pendulum on an 
n — 1-dimensional sphere S""^ as an example of a mechanical second-class constraints system and 
the 0(n) non-linear sigma model as an example of a field theory under second-class constraints are 
discussed in details and quantized using the existence of underlying dilatation gauge symmetry and 
by solving the constraint equations explicitly. The underlying gauge symmetries involve, in general, 
velocity dependent gauge transformations and new auxiliary variables in extended configuration 
space. Systems under second-class holonomic constraints have gauge-invariant counterparts within 
original configuration and phase spaces. The Dirac's supplementary conditions for wave functions of 
first-class constraints systems are formulated in terms of the Wigner functions which admit, as we 
show, a broad set of physically equivalent supplementary conditions. Their concrete form depends 
on the manner the Wigner functions are extrapolated from the constraint submanifolds into the 
whole phase space. 

PACS numbers; ll.lO.Ef, ll.lO.Lm, 11.15.-q, 11.30.-j, 12.39.Fe 



I. INTRODUCTION 

The specific feature of gauge theories is the occurrence 
of constraints which restrict the phase space of the system 
to a submanifold. A systematic study of the Hamiltonian 
formulation of gauge theories was made by Dirac 0, Q 
who classified the constraints and developed the opera- 
tor quantization schemes of the constraint Hamiltonian 
systems. 

The Dirac's theory of constraint systems was com- 
bined further with the Feynman path integral method. 
Faddeev jSj found an exphcit measure on the constraint 
submanifolds, entering the path integrals. The Feyn- 
man diagram technique for non-abelian gauge theories 
was developed by Faddeev and Popov Q. The effective 
QCD Lagrangian involving ghost fields obeys the Becchi- 
Rouet-Stora (BRS) symmetry 0|. A completely general 
approach to quantization of gauge theories, in which the 
BRS transformation acquires an intrinsic meaning, is de- 
veloped by Fradkin and his collaborators jfl,!?!. An alter- 
native symplectic scheme has been proposed by Faddeev 
and Jackiw |^ . An introduction to quantization of the 
gauge theories with the use of the path integral method 
can be found in Refs. [13, 113. The path integral repre- 
sentation for the evolution operator satisfies the unitarity 
condition and meets requirements of the relativistic co- 
variance. 

According to the Dirac's classification, constraint 
equations like Q.a = 0, appearing in the gauge theories, 
are of the first class. The Poisson bracket of first-class 
constraint functions is a linear combination of constraint 
functions 

{nA,nB]^CA^nD (LI) 



where Cj^g are some structure functions. The quantiza- 
tion of gauge theories involves a set of gauge- fixing condi- 
tion xa = 0. The functions XA must be admissible, i.e., 
the Poisson bracket of the gauge-fixing and constraint 
functions is non-degenerate: 

det{xA,ris} ^0. (1.2) 

The gauge-fixing functions are independent 

{xa,Xb} = Q. (1.3) 

Physical observables are gauge invariant and do not de- 
pend on the choice of xa- 

The wave functions satisfy the operator Dirac's sup- 
plementary condition 

(Ib"^ = 0. (L4) 

Constraint functions Qa of second-class constraints sys- 
tems obey the Poisson bracket relations which form a 
non-degenerate matrix 

det{g„ 06)^0. (L5) 

The number of the constraints is even (a = 1, 2m). 

Among physical systems which belong to the second- 
class constraints systems are anomalous gauge theories 
m mill m Hg, the 0(4) non-Hnear sigma model, 
which constitutes the lowest order of the chiral perturba- 
tion theory [ITi] , many-body systems involving collective 
and independent-particle degrees of freedom 18j, and 
others. 

The Dirac's quantization method of such systems con- 
sists in constructing operators reproducing the Dirac 
bracket for canonical variables and taking constraints to 



2 



be operator equations of the corresponding quantum the- 
ory. 

Batahn and Fradkin [T^ developed a quantization pro- 
cedure for the second-class constraints systems, which 
converts constraints to the first class by introducing new 
canonical variables. The problem reduces thereby to the 
quantization of a first-class constraints system in an en- 
larged phase space. This method (see also ^J) was found 
to be particularly useful for construction of the effective 
covariant Lagrangians in an extended configuration space 

mill. 

The Hamilton- Jacobi scheme is also used for the quan- 
tization of constraint systems [2II2II2I I2ll2^ . 

Gauge systems are quantized by imposing gauge-fixing 
conditions. The initial system is reduced to a second- 
class constraints system. The evolution operator in the 
path representation can be written as 0, 0| 

/in, -JTL 

xexpl^^Jdt{fq,-H)^ (L6) 

where Ti. is the Hamiltonian. The path integral represen- 
tation ljl.6l) allows not to solve the constraint equations 
explicitly and work in the unconstrained phase space. 

The usefulness of the reduction of second-class con- 
straints systems to equivalent gauged systems consists in 
getting the supplementary condition (|I.4|I which is not 
evident otherwise. 

The second-class constraints systems to which gauge 
systems are reduced obey specific requirements: 

(i) The constraint functions Ga split naturally into 
canonically conjugate pairs Qa ={xa, ^b), 

(ii) the wave functions satisfy Ea. (|I.4p . 

(iii) the gauge-fixing functions xa are identically in 
involution Ea. (|I.2|) . and 

(iv) the constraint functions CIb associated to the 
gauge invariance are first class Ea. (|l.l() . 

As a consequence of Eg. ljLSI) . the matrix {xa,^b} is 
non-degenerate. The constraint functions are not defined 
uniquely. In particular, a linear transformations of Qa by 
a non-degenerate matrix leads to constraint functions Q'^ 
describing the same constraint submanifold. 

A question arises if constraints of an arbitrary second- 
class constraints system can be redefined to fulfill (i) - 
(iv)? 

This is the case for holonomic systems 27]. Such sys- 
tems can be treated as being obtained upon a gauge- 
fixing of the corresponding gauge invariant systems. 
Within the generalized Hamiltonian framework, con- 
straints of holonomic systems are of the second class. In 
the Lagrangian formalism, the corresponding constraints 
do not depend on generalized velocities. 

In this work, we review the underlying gauge symme- 
tries of the holonomic systems and report new results 
connected to the quantization of more general second- 
class constraints systems. 



The paper is organized as follows: In the next Sec- 
tion, we start from discussing a simple, but instructive 
example of a mechanical system under second-class holo- 
nomic constraints. The one-dimensional reduction of the 
0{n) non- linear sigma model is discussed, which is equiv- 
alent to a mathematical pendulum on n — 1-dimensional 
sphere in an n-dimensional Euclidean space. Lagrangian 
C of the system is transformed on the constraint sub- 
manifold of the configuration space into an equivalent 
Lagrangian to make explicit the appearance of an un- 
derlying gauge symmetry. The corresponding Hamilto- 
nian, its constraints structure, and transformation prop- 
erties are described in Sect. 3. 

This example is analyzed further in Sect. 6 to con- 
struct a link with the Dirac's quantization method. 

In Sects. 4, 5, and 8, the quantization of systems under 
the second-class constraints is discussed. 

The Poison bracket of the constraint functions Qa 
forms a symplectic structure in the space of constraint 
functions. The corresponding symplectic basis is suit- 
able for splitting the constraints into canonically conju- 
gate pairs {xa,^a)- In Sect. 4, an algorithm is proposed 
for constructing the functions XA and U,a, describing the 
constraint submanifold of second-class constraints sys- 
tems, for which the involution relations {xAtXb} = 0, 
{Q.a,^b} = 0, and {xa,^b} = 5ab hold in the strong 
sense in an entire neighborhood of any fixed point of the 
constraint submanifold. 

In Sect. 5, properties of the gradient projections into 
constraint submanifolds described by equations xa = 0, 
Q.A = 0, and onto their intersection are discussed. The 
gradient projection is useful for calculation of the Dirac 
bracket and for constructing symplectic basis for the con- 
straint functions. 

The Dirac's quantization of the model of Sect. 2 is 
performed in Sect. 6 and the equivalence with the results 
of Sect. 3 is established. In Appendix A, the model of 
a mathematical pendulum on n — 1-dimensional sphere 
is embedded further into a space of higher dimension to 
study the constraints structure and the associated gauge 
invariance. The equivalence with the results of Sects. 3 
and 6 is proved. 

In quantum mechanics, systems are described by wave 
functions or, equivalently, by density matrices or the 
Wigner functions. The Wigner functions are the most 
closely related to the classical probability densities in the 
phase spaces. In Sect. 8, we analyze constraints to be im- 
posed on the Wigner functions within the framework of 
the generalized Hamiltonian dynamics. Such constraints 
are not unique and depend decisively on the way the 
Wigner functions are extrapolated from the constraint 
submanifold into the unconstrained phase space. The 
main result of this Section consists in demonstrating the 
fact that constraints imposed on the Wigner functions 
can be taken in a symmetric fashion with regard to per- 
mutations XA ^ ^A, XA ^ —^A and, more generally, 
canonical transformations preserving the Poisson brack- 
ets. Those constraints are not symmetric in the Dirac's 
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quantization scheme. We show that transition ampH- 
tudes are not affected by canonical transformations mix- 
ing the constraint functions. 

In Sect. 9, the 0(4) non-hncar sigma model is dis- 
cussed as a field theory counterpart of the mathemati- 
cal pendulum. We perform quantization in the straight- 
forward way by solving the constraints from the out- 
set and demonstrate the equivalence of the results with 
the method based on construction of the gauge-invariant 
model of the 0{4) non-linear sigma model. The La- 
grange measure for the 0(4) non-linear sigma model, 
entering the path integral in the configuration space, is 
constructed and a parameterization is proposed for the 
pion fields, in terms of which the perturbation theory 
consistent with the mean field (MF) approximation can 
be developed. The covariance and the unitarity of the 
^-matrix are demonstrated. 

The results are summarized in Sect. 10. 



II. GAUGED LAGRANGIAN FOR A 
SPHERICAL PENDULUM 

We start from discussing a simple example of a me- 
chanical system under holonomic constraints: A spher- 
ical pendulum on an n — 1-dimensional sphere S'^~^ in 
n-dimensional Euclidean space. The trajectories of the 
system are determined as conditional extremals of the ac- 
tion functional A = J Cdt on the constraint submanifold 
X = 0, with 

X = ln(j) (II.l) 
where (fy^ = and a = 1, n. Lagrangian 

c = (n.2) 

together with the constraint x — defines a pointlike 
massive particle on an n — 1-dimensional sphere S"^^ 
as a mechanical analogue of the 0{n) non- linear sigma 
model. The constraint x = implies 

= 0, (II.3) 
+ 0"^" = 0, (II.4) 

The dots stand for higher time derivatives of the con- 
straint function x- Ea. (jII.3p shows that the radial compo- 
nent of the velocity vanishes. The second equation shows 
that the radial component of the acceleration equals to 
the centrifugal force. 

The equations of motion can be found using the 
D'Alembert-Lagrange variational principle for condi- 
tional extremals of the action functionals, or equiva- 
lently, the Euler-Lagrange variational principle for un- 
conditional extremals with the constraints implemented 
through the Lagrange's multipliers method. 

Let us substitute £ — > £i = C+Xx and perform uncon- 
ditional variations over the Lagrange's multiplier A and 



the coordinates (j)"- One gets x = and 0" = Xcj)" / (j)'^ . 
Multiplying further the last equation by and substi- 
tuting the result into Ea. ljIL4l) . one gets A — — 
Given that the A is fixed from Ea. l|IL4|l . the higher time 
derivatives of x vanish identically. The radial component 
of the acceleration is determined by the constraint, while 
the tangent component of the acceleration vanishes. The 
equations of motion can be presented in the form 

A'^^{cf>)^icj^t'/cj^)=0 (n.5) 

where 

A"'3(0) = S"^ - (t>°'(l)'^ /(j)'^. (IL6) 

The tensor defined above obeys 

0"A"^(</)) = 0, (n.7) 
A"'3(<^)A/57(0) ^ A"^(0). (n.8) 

It is invariant also with respect to dilatation 

A"'5((^') = A"'3(0), (II.9) 

for 

0" ^ 0'" = exp(6i)</)" (11.10) 

where 9 is an arbitrary parameter. 

Eas. (jII.5|) tell that the particle moves without tangent 
acceleration. In general, the acceleration orthogonal to 
a constraint submanifold x = is fixed to keep the par- 
ticle on it. In our case, the radial component of the 
acceleration is determined by the constraint. Eas. (|II.5|) 
can also be derived using the D'Alembert-Lagrange vari- 
ational principle. 

On the submanifold x = 0, the radial components of 
the velocities are equal to zero Ea. ljII.3|) . so one can re- 
place (/)" by the tangent velocities A"''((/))0^ in C. 

The conditional extremals of the action functional A = 
J Cdt do not change also if we divide Chy cf)^ {— 1). The 
conditional variational problem for Lagrangian 

= iA"/'(0)r(/./'/02 (11^^^) 

is thus completely equivalent to the conditional varia- 
tional problem we started with. 

The equations of motion (|II.5p determine uncondi- 
tional extremals of the action functional A^: = ^ C^^dt 
on the configuration space M ~ {4>")- 

The extremals of the action functionals A and ,4* un- 
der the constraint x = coincide. 

At, depends on the spherical coordinates (p^/c/) which 
lie on an n — 1-dimensional sphere of a unit radius. 
Eas. ljII.5|l for unconditional extremals of the A^, in the 
coordinate space cj)" coincide, as it should, with the 
D'Alembert-Lagrange equations for extremals of A^, in 
the space formed by the spherical coordinates under the 
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condition that they belong to an n — 1-dimensional sphere 
of a unit radius. 

It is seen that Lagrangian is invariant with re- 

spect to dilatation Ijll.lOp where 6 is an arbitrary function 
of time. In the context of the dynamics defined by HII.11() 
with no constraints imposed, (j) turns out to be an arbi- 
trary function of time. It can always be selected to fulfill 
the constraint x = or some other admissible constraint. 

The constraint X — ^ can therefore be treated as a 
gauge-fixing condition, the function (f> as a gauge degree 
of freedom, the ratios 4>°'/4> o,s gauge invariant observ- 
ables. 

The equations of motion Ijll.Sp are formulated in terms 
of the gauge invariant observables. 

The above gauge symmetry is defined "on-sheU", i.e., 
for (/)" treated as a time derivative of 0". In the tangent 
bundle TA^ = where the coordinates and their 

derivatives are independent, is invariant with respect 
to a two-parameter set of transformations: 



and 



6'" = exp(6')(/.", 
6'" = exp(6')(/." 

r = r, 

l/a la. I ^ j,Q 



(11.12) 
(11.13) 



(11.14) 
(11.15) 



The last two equations describe the invariance under vari- 
ation of the radial component of the particle velocities. 
If the 0" is treated on-shell, then e ^ 6. 

A pointlike particle on an n — 1-dimensional sphere 
S"^^ has therefore underlying gauge symmetry con- 
nected with dilatation of the coordinates 0". Its physi- 
cal origin is simple: We allow virtual displacements from 
the constraint submanifold and treat such displacements 
as unphysical gauge degrees of freedom. The constraint 
equation x = turns thereby into a gauge-fixing condi- 
tion. The physical variables are specified by projections 
of the coordinates (fi" onto an n — 1-dimensional sphere 
of a unit radius. Those projections are the spherical co- 
ordinates (jf /4>- 

Let us consider the system (jll.ll|) within the general- 
ized Hamiltonian framework. 



III. 



THE GAUGED SPHERICAL PENDULUM 
HAMILTONIAN 



Lagrangian Ijll.lip is defined outside of the constraint 
submanifold x — ^ ^^^id invariant with respect to the 
dilatation (|II.10|I . The value (f) is an arbitrary function 
and becomes a gauge degree of freedom. The equations 
of motion (|II.5|I are derived using the constraint x = 
0. Being formulated, the equations of motion do not 
depend, however, on the constraint anymore and allow 
an extension to the unconstrained configuration space 
Ai. The invariance under the dilatation, 



consequence of the two-parameter set of global symmetry 
transformations, Eas. (|11.12l) - (|11.15() . of as a function 
defined on the tangent bundle TAA. 

One can consider therefore A<t without imposing any 
constraints and treat the equation x = as a gauge-fixing 
condition. 



A. Gauged spherical pendulum within generalized 
Hamiltonian framework 



The canonical momenta corresponding to (jf are de- 
fined by 

7r" = ^ = A"'5(0)//</)2. (III.1) 

They satisfy constraints 

tt" - A"'3(0)7r'3 «0 (IIL2) 

which are equivalent to the one primary constraint: 

= 07r « 0. (III.3) 

The primary Hamiltonian can be obtained with the use 
of the Legendre transform: 



(IIL4) 



For n = 3, 7i is proportional to the orbital momentum 
squared. The non-vanishing Poisson bracket relations for 
the canonical variables are defined by 



(IIL5) 



The constraint function is stable with respect to the 
time evolution: 



The relations 



{n,n} = Q. 



(III.6) 

(IIL7) 
(IIL8) 



show that n generates dilatation of (p" and tt". 
The transformation law for the canonical coordinates 
Ea. (|III.7p is in agreement with Ea. (|II.12|) in its infinites- 
imal form. The transformation law for the canonical 
momenta, which follows from Eas. l|II.12p . ljII.13|) . and 

(dnj, 

tt" ^ tt'" = cxp(-6l)7r", (IIL9) 

considered in its infinitesimal form, is in agreement with 
Ea. (|III.8p either. The Hamiltonian Ti is gauge invariant 
under the dilatation. 

The roles of the gauge-fixing function x and the gauge 
generator are similar. The function x is identically in 
involution with the Hamiltonian: 



{x,n} = o. 



(IH.IO) 
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The Poisson bracket relations 

= 0, (III.ll) 

{^",x} = (ni.i2) 

define the one-parameter set of transformations with re- 
spect to which Ti is invariant. The function x generates 
shifts of the longitudinal component of the canonical mo- 
menta. This symmetry is connected to the invariance of 
described by Eas. HII.14| and (jll.lSp . 
The gauge-fixing condition x = is admissible: 

{xM = ^- (111.13) 

The equations of motion generated by the primary 
Hamiltonian look like 

tt" {7r",7^} = -0"A''T(<?!))7r^7r'^. (III.15) 



B. Quantization of spherical pendulum 

The quantization of a mathematical pendulum on an 
n— 1-dimensional sphere is discussed in Ref. [22] where the 
standard Batalin-Fradkin-Tyutin (BFT) algorithm 0, 
for second-class constraints systems is applied. The 
second-class constraints appear if one starts directly from 
C and formulate the conditional variational problem for 
X = 0. By constructing auxiliary fields, it is possible to 
pass over to an equivalent first-class constraint system. 

In our approach, we start from As^ which is gauge in- 
variant explicitly, so the constraints appear to be of the 
first class from the start. One can therefore quantize the 
pendulum as a gauge-invariant system without introduc- 
ing auxiliary fields. 

From the point of view of the generalized Hamiltonian 
framework, the gauge-fixing conditions and the gauge 
generators play similar roles. The function x does, how- 
ever, not generate transformations in TM.^ so it appears 
just as a candidate for gauge-fixing function. If we pass 
to the Lagrangian framework, we can verify that x = 
is the gauge-fixing condition indeed. 

The system has the gauge invariance described by the 
generator and the admissible gauge-fixing condition 

x = o. 

The standard procedure for gauge theories can therefore 
he applied for quantization of the spherical pendulum. 

The system is quantized by the algebra mapping 
((/)", tt") ^ (0",7r") and {,} ~> -i/hl]. Consequently, 
to any symmetrized function in the phase space variables 
one may associate an operator function. The function is 
symmetrized in such a way that quantal image is a her- 
mitian operator. 

The quantum hermitian Hamiltonian has the form 

n = i(/)A°''7r'3(^A"^7r^. (IH.16) 



The vector icjjA^^n^ gives the angular part of the gra- 
dient operator. Although it is not conspicuous, H does 
not depend on the radial coordinate 0. The constraint 
operator can be defined as 

n = (0"7r" + TT°'(t>")/2. (IH.17) 

It acts only on the radial component of (j)" , so the relation 

[n,n]=0 (IH.18) 

holds. 

The physical subspace of the Hilbert space is singled 
out by imposing the Dirac's supplementary condition 

= 0. (III.19) 

This condition implies 

* = 0-"/2*i(0"/0). (IH.20) 

The physical information is contained in / (j)) . 

The path integral for the evolution operator becomes 

z=j n(£|^'^(^)^(")-p{^ / 
(1H.21) 

Eas. (|III.19l) and ljIII.2ip solve the quantization prob- 
lem for a mathematical pendulum on an n — 1- 
dimensional sphere 5""^. 

It is desirable that the quantization procedure does 
not destroy the classical symmetries which results in hav- 
ing the supplementary condition (|III.19p satisfied by the 
state '^{t) for any value of t. This feature can, in gen- 
eral, be violated either due to complex terms entering the 
Hamiltonian or by approximations adopted for treating 
the operator eigenvalues. 

Given the initial state wave function ^(0) satisfying 
(|III.19ll . the final wave function '^{t) can be found ap- 
plying the evolution operator ljIII.21|l on ^"(0). Since 
VL commute with the Hamiltonian, the final state wave 
function obeys Eas. HIII.19|l . 

In the standard canonical frame 0|, the first canoni- 
cal coordinates and momenta are identical with the con- 
straint functions. If the Poisson bracket of the Hamilto- 
nian with the constraint functions vanishes, the Hamilto- 
nian is independent of the first canonical coordinates and 
momenta. It means that the quantized Hamiltonian is 
commutative with operators associated to the constraint 
functions. 

By initiating the quantization procedure in the stan- 
dard canonical frame, one gets the classical symmetries 
preserved and the Dirac's supplementary condition ful- 
filled at any time. 

The fact that the dilatation symmetry for the spherical 
pendulum is preserved on the quantum level is connected 
with the fact that the canonical frame we used is simply 
related to the standard canonical frame. 
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The integral over the canonical momenta Ea. (|III.21|) 
can be simplified to give 



where 



Z 



n 



(9x/5r)''5(x)rf>exp<{^ 



(III.22) 

Lagrange's measure \l {dx/ d(j)°')^ 5{x)d''^4' coincides with 
the volume element of S*"^^ sphere. It can be rewritten 
as an invariant volume of the configuration space, e.g., 
in terms of the angular variables ((/?!. with the 
help of the induced metric tensor. It is invariant under 
0(n) rotations and remains the same for all functions 
X vanishing aX <j) — 1. For n = 4, Lagrange's measure 
matches Haar's measure of the group SU{2). 

In gauge theories, the evolution operator is indepen- 
dent of the gauge-fixing conditions [Sj. We can insert 
det{x, f^} = 1 into the integrand of Ea. HIII.22() to bring 
the measure into the form identical with gauge theories. 
In the path integral, x can then be replaced with an arbi- 
trary function. The condition ljIII.19|l does not comprise 
the constraint x = also. 

It is remarkable that physical observables depend on 
half of the number of the second-class constraints. 

The main question to be raised here if the quantization 
method described above is general enough or specific only 
for the spherical pendulum? In Sect. 7 we show that such 
a method works for mechanical systems and field theories 
under holonomic constraints. The quantization of more 
general second-class constraints systems is discussed in 
Sect. 8. 



IV. LOCAL SYMPLECTIC BASIS FOR 
SECOND-CLASS CONSTRAINTS FUNCTIONS 

Two sets of the constraint functions are equivalent 
if they describe the same constraint submanifold. The 
Hamiltonian function admits transformations which do 
not change its value and its first derivatives on the con- 
straint submanifold. This allows to make transforma- 
tions of the constraint and Hamiltonian functions with- 
out changing the physical content of theory. 

The second-class constraints satisfy det{Qa,Gb} 7^ 0. 
The Poisson bracket defines therefore a non-degenerate 
symplectic linear structure in the vector space of the con- 
straint functions Ga- Indeed, any linear transformation 
Q'g^ = OabGb with matrix Oab depending on the canoni- 
cal variables transforms accordingly the Poisson bracket: 
{Ga,Gt} « OacObd{Gc,Gd}- The Poisson bracket plays 
thereby a role of a skew-scalar product in the symplectic 
vector space of the constraint functions. Every symplec- 
tic space has a symplectic basis (see, e.g., |23|), so the 
constraint functions can be brought by linear transfor- 
mations into the form 





— En 



E„ 




(IV.2) 



with Em being the mxm identity matrix, Xablbc — ~Sac 
Using representation Ga = {xa,^a), one has 



{flA^nB} « 0, 
{xa,Xb} « 0, 
{xa,^b} « Sab- 



(IV.3) 
(IV.4) 
(IV.5) 



This basis is not unique. Indeed, there remains a group of 
symplectic transformations Sp{2m), which keeps the Poi- 
son bracket of the constraint functions in the symplectic 
form. 

At any given point of a neighborhood of the constraint 
submanifold one can find symplectic basis for second-class 
constraints functions in the weak form. 

This result can be strengthened as shown below: 

(„4) If ^ is close enough to the constrain submani- 
fold, one has det{xA{0:^B{£,)} 0- By continuity 
there exists a finite neighborhood of such that 
det{xA{C),^B{C)} 7^ V^' G Aj. Let us assume 
that the intersection of A^ with the constraint subman- 
ifold is not empty, i.e., one can find ^1 G A^ such that 
Ga{£,i) = 0. If it is not fulfilled, we start from another ^ 
closer to the constraint submanifold. 

Let us chose an equivalent set of the constraint func- 
tions XA ^ x'a = {xa,^b}~^Xb to ensure {x'a,^b} « 
Sab in the region A^. 

{B) We replace further Qa ^'a = ^A — \Cabx'b to 
get equations 

{^'aMb} ~ {^aM-Cab + \cadCbfWd.^f} « 0. 

(IV.6) 

These equations can be solved for matrix Cab in terms 
of a power series of the matrix {x'd, x'fY- 



where C^ab 



Cai 



-C^BA 



fc=0 



(IV. 7) 



C 



[0] 
AB 



{^A.^b}, 



C 



[k+l] _ 1 



AB 



4 ^ 

J=0 



^AD '^BFIXdtXfU 



with fc = 0, 1, ... . 

(C) The transform x'a x'a ^ {x'a^^'b}~^x'b brings 
back the Poisson bracket {x'a, ^'b} to the diagonal form 
Sab- 

iV) The last transform looks like x'a ~^ x'a — x'a ~ 
Hx%XbWb- It provides {xa^x'b} « Oand keeps the 
relation {x'a,^'b} ~ Sab unchanged. 

Now, wc remove primes from the notations. As a result 
of the steps (y^)-(P), we obtain weak equations 



{Ga, Gb} ~ lab 



(IV.l) 



{GaiO^QbiO} 



(IV.8) 
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VC S A^. It is manifest that Ea. (|IV.8p is valid in some 
neighborhood of any point of the constraint submani- 
fold too. 

The symplectic basis for second-class constraints func- 
tions in the weak form exists in an entire neighborhood 
of any given point of the constraint submanifold. 

The existence of the local symplectic basis in the weak 
form is on the line with the Darboux's theorem (see, e.g., 
) which states that around every point ^ in a symplec- 
tic space there exists a coordinate system in such that 
^ S Aj where the symplectic structure takes the standard 
canonical form. The symplectic space can be covered by 
such coordinate systems. 

This is in sharp contrast to the situation in Ricman- 
nian geometry where the metric at any given point x 
can always be made Minkowskian, but in any neighbor- 
hood of X the variance of the Riemannian metric with 
the Minkowskian metric is, in general, ~ Ax^. In other 
words, by passing to an inertial coordinate frame one can 
remove gravitation fields at any given point, but not in 
an entire neighborhood of that point. The Darboux's 
theorem states, reversely, that the symplectic structure 
can be made to take the standard canonical form in an 
entire neighborhood Aj of any point ^ € Aj. In Rie- 
mannian spaces, locally means at some given point. In 
symplectic spaces, locally means at some given point and 
in an entire neighborhood of that point. 

Locally, all symplectic spaces are indistinguishable. 
Any submanifold in a symplectic space, including any 
constraint submanifold, is a plane. The possibility of 
finding the standard canonical frame [s^l illustrates this 
circumstance. 

In the view of this marked dissimilarity, the validity of 
Eas. ljIV.8p in a finite domain looks indispensable. 

The global symplectic basis exists apparently for m = 
1 as, e.g., in the case of the spherical pendulum. The 
global existence of the basis (|IV.8|I has been proved for 
systems with one primary constraint |29|. and also, as- 
suming that det{xA,^B} 7^ holds globally |30| . 

Admissible transformations on the second-class con- 
straints functions allow to bring Eas. (|IV.8|) into a strong 
form. The Hamiltonian can also be modified to convert 
the Poisson bracket relations with the second-class con- 
straints functions into the strong form without changing 
the dynamics. The arguments given below follow closely 
our discussion of holonomic systems [27j : 

According to Eas. ljIV.8|) . at any given point of the con- 
straint submanifold where det{Oa, Qb} 7^ one can select 
symplectic basis in which the constraint functions satisfy 

{ga,gb}=Iab+C,'bGc (IV.9) 

in an entire neighborhood of that point. The first-class 
Hamiltonian Ti. has relations 

{Ga^nj^n^Gb- (iv.io) 

The Jacobi identities for Qa, Gb, and Gc and for Ga, Gb, 



and Ti. imply 

+ Cbca+Ccab « 0, (IV. 11) 

TZab-Uba « (IV. 12) 

where Cabc = 1c.dC J and Tlab = laJla ■ 
Let us define 

G'a = Ga + \c^'GbGc, (IV.13) 

W = n-]^n''^GaGb- (IV.14) 

The constraint functions Ga and G'a coincide on and in 
vicinity of the constraint submanifold up to the second 
order in Ga and generate accordingly identical phase flows 
on the submanifold Ga = ^- The Hamiltonian functions 
TL and TC coincide on and in vicinity of the constraint 
submanifold up to the second order in Ga and generate 
on the submanifold Ga = ^ identical Hamiltonian phase 
flows. 

Using Eas. pV.Up and l|IV.12p one gets 

{G'a,G'b} = lab + C^iKG'a. (IV. 15) 

{G'a.n'} = n^'^G'bG', (IV.16) 

where Cjf' and T^-j*^ are new structure functions. The 
first-order terms in Ga do not appear in the right sides 
of these equations. The pair (^^,W) describes the same 
Hamiltonian dynamics as (GajTi), being at the same time 
in a stronger involution around the constraint submani- 
fold. 

The above procedure can be repeated to remove the 
quadratic terms in Gai cubic terms in Ga, etc. In general, 
assuming 

{GI^KgI'^} - Iab + Cj^-'''Gif...Gi':}, (IV.17) 
= TZa'^-'^Gi'^-Gl!'^, (IV.18) 

we get as a consequence of the Jacobi identities and of 
the symmetry of structure functions C^^^ "'^'' and TZ^^ - -^'' 
with respect to the upper indices: 

Cabci...Ck ~^ Cbcia...Ck ^cia6...Cfc ~ 0? (IV. 19) 

nab^...b,~nb,a...b, « 0. (IV.20) 

The next-order constraint functions and the Hamiltonian 
are given by 

^[fc+1] ^ ^[fc] —n^0bl...b,glk]g[k] glk]^^^y.^^^ 

k ~\~ 1 1 fc 

Using Eas. ljIV.17p and (jIV.18j) . one can calculate the 
next-order structure functions ^nd TZa^ " ''^^ 

and repeat the procedure. If structure functions vanish, 
we shift k by one unit and check Eas. (jIV.17j) and ljIV.18j) 
again. At each step, the Poisson bracket relations get 
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closer to the normal form. In the limit k — > +00, we 
obtain 

{Ga.Gb} = lab. (IV.23) 

{Ga.n} = (IV.24) 

where Ga = limfe_oo Sl''' and H = limfc^oo H^'"!. The 
matrix Jab defines splitting of the constraints into two 
groups (xaj^a), such that 

{xa,Xb} = 0, (IV.25) 
{hA.hB} = 0, (IV.26) 
{xaM = 5ab. (IV.27) 

The progress made consists in extending the valid- 
ity of Eas. l|IV.ip from one point into its neighborhood 
Eas. ljIV.8|l . and further, in passing from the weak to the 
strong form of the Poisson bracket relations Eas. ljIV.23p . 

The symplectic basis for second-class constraints func- 
tions exists in the strong form in an entire neighbor- 
hood of any given point of the constraint submanifold. 
There exits a Hamiltonian function TL, describing the 
same dynamics on the constraint submanifold as the ini- 
tial Hamiltonian function H, which is identically in in- 
volution with the constraint functions. 

An independent geometric-based construction of Ga 
and Ti. is given in Sect. 5. 

A similar local basis exists for first-class constraints 
0,|33- The arguments of Refs.0,|33 apply to second- 
class constraints under special restrictions which are dis- 
cussed in Appendix B. 

For systems of pointlike particles under holonomic con- 
strains, Eas. HIV.9p hold globally, so Ga and H. exist glob- 
ally also. Furthermore, VLa is linear in the canonical mo- 
menta, xb does not depend on the canonical momenta, 
and Ti. splits into a sum of a kinetic energy term quadratic 
in the canonical momenta and a potential energy term 
depending on the canonical coordinates 0| . 



V. GRADIENT PROJECTION 

The concept of the gradient projection is useful for 
applications. It defines functions ^s(C) which project an 
arbitrary point ^ of the phase space onto a submanifold 
= (a = l,...,2m) of the phase space along phase 
flows generated by the constraint functions Ga- 



A. Full gradient projection 

Let = be second-class constraints, det{Ga,Gb} 7^ 
0, ^* {i = l,...,2n) are canonical variables. In vicinity 
of the submanifold Ga = 0, the projections can be con- 
structed explicitly. Near the constraint submanifold, one 
can write = ? + {^iGa^^a- The small parameters 



Aq are determined by requiring Ga{S,s{C)) = to the first 
order in Ga- We get 

UO=i-{i,Ga}\\{Ga,Gb}\\-'Gb- (V.l) 

It is seen that {^?a,Cs(0} = consequently 
{^a,/(Cs(0)} = fo'^ any function /. This is natural, 
since Ga generate phase flows along which the projections 
^s(^) have been constructed. 

The reciprocal statement is also true: If {^a,/} — 
for all Ga, then / = f{£,s{£,))- Indeed, the coordinates 
on the constraint submanifold can be parameterized by 

- The coordinates describing shifts from the constraint 
submanifold can be parameterized by Ga- The functions 
/ can in general be written as f = f{^s,Ga)- If the 
Poisson brackets of / with all Ga vanish, / depends on 
only. 

This can be summarized by 

{GaJ}^0^f = f{UO)- (V.2) 

Beyond the lowest order in Ga , the operation is unique 
provided the phase flows commute. This is always the 
case for the constraint functions taken to accomplish 
(|IV.23|) in a finite neighborhood of the constraint sub- 
manifold. 



B. Partial gradient projection 

Let Ga split into XA and Qa- We wish to construct 
functions which project an arbitrary point ^ of the 

phase space onto the gauge-fixing surface XA = with 
the use of the constraint functions Qa associated to gauge 
transformations. In vicinity of the submanifold, one can 
write £,u{£,) = C + ^a}^a- To the first order in XA, the 
parameters can be found from equation Xa{£,u{0) — 
0: 

UO = (.-{t^A}\\{xA,nB}\\-\B- (V.3) 

The projection is made along the phase flows of Qa, 
so {Ha, f{£.u{£.))} = for any function /. 

The reverse statement is also true. We write / = 
f{^u,XA) and conclude from {17^,/} = that the de- 
pendence on XA drops out. 

It can be summarized as follows: 

{nAj}^o^f^f{UO)- (V.4) 

The second partial gradient projection can be made 
onto the submanifold flA = with the use of the con- 
straint functions XA- The result is similar to Ea. l|V.3p 

UO^^ + {tXA}\\{XAMB}\\-'nB- (V.5) 

The relation {xA, /(Cti(O)} = is valid for any func- 
tion /. Furthermore, 

{XA,/} = 0^/ = /(e„(O). (V.6) 
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Combining the partial projections, e.g., £,s{0 = 
iv{£,u{0)j 0116 gets the full gradient projection. To the 
first order in Qa, the order in which the partial projections 
are apphed does not matter, so £,v{^u{£,j) = Cu{£,v{0)- 
The full gradient projection constructed in this way co- 
incides with that given by Ea. (|\rT)l . 



C. Example: Gauged spherical pendulum 

The constraint function f2 can be used to bring the 
vector onto a sphere of a unit radius 4> — I (x ^ 
0). Such a transformation has a meaning of a gradient 
projection. The functions (j)" and tt" can be constructed 
in terms of the variables 0" and iTa- 

C = cxpi9)r, (V.7) 
< = cxp(-6l)7r". (V.8) 

The condition xi't'in'^u) = gives exp(6') — 

The Poisson bracket relations for the projected vari- 
ables (|V.7p and HV.8|) can be found to be 

{ruAi} = 0, (V.9) 

{<^S,7rf} = K'{cl>u). (V.IO) 

= €<-<l>Z^i (v.ii) 

where A^'^ = 5°^P - K<l>t 

The following properties are worthy of mention: 

(i) The Poisson bracket of the projected canonical vari- 
ables coincides with the Dirac bracket associated to the 
constraints x = ^-nd 51 = on the submanifold x = 0- 

The Dirac bracket for the canonical variables can be 
calculated using Ea. HV.25p to give 

{r,^^}D = 0, (V.12) 
{r,^^}D = A^^{<i>), (V.13) 

{t:",t:'^}d = ((/)'^^"-(/)"^'^)/(/)'. (V.14) 
The right sides of Eas. ljV.9p - (|V.11|) are reproduced at 

x = o. 

(ii) The relations IIV.9|) - IIV.ll|l define a Poisson al- 
gebra in the space of functions /(0",7i'") depending on 
the projected canonical variables, so they can be used 
to generate consistently a Hamiltonian phase flow on the 
constraint submanifold x = 0. 

(iii) The Hamiltonian fmiction 

H=^AZP{cl>u)K^t (V.15) 

coincides with Eq JTTLlll : H= n{(j)Z,K)^ H(0",7r"). 
The Hamiltonian ljHI.4p is thus the function of the pro- 
jected variables S,u{0- I* can be defined first on the sub- 
manifold X = s-nd then extended to the unconstrained 
phase space using the gradient projection parallel to the 
phase flow associated to 57. The relation {57, 7i} = 0, 
Ea. (|III.6|l . is the necessary and sufficient condition (see 



Ea. (jV.4|) ') for 7i to be a function of a fewer number of 
variables H = 

Let us consider the gradient projection onto the sub- 
manifold f2 = using the constraint function x- The 
constraint x — is responsible for shifts of the longitu- 
dinal component of the canonical momenta Eas. (|lll.ll() 
and piI.12|) . The functions 0" and tt" have the form 

C = (V.16) 
< = tt" -0"<?!)7r/(/.2. (V.17) 

Equation 57 ((/>„, tt^) = is fulfilled identically. 

The Poisson bracket relations for the projected vari- 
ables (|V.16|) and (|V.17|) can be found to be 

{C'^i^} = 0, (V.18) 
{r.y,} = A'^^i^,), (V.19) 
{<,7r,^} = (<^e<-0>.^)/0'- (V.20) 

One can see again: 

(i) The Poisson bracket of the projected variables co- 
incides with the Dirac bracket Eas. ljV.12|) - ljV.14|) as- 
sociated to the constraints x = ^'^d 57 = on the 
submanifold 57 = 0. 

(ii) The Poisson bracket relations (jV.lSp - ljV.20p are 
closed and define thereby a Poisson algebra in the space of 
functions depending on the projected canonical variables 

(iii) The Hamiltonian function 

n = i0?,A"''(0.)<^^ (V.21) 

coincides with Eq.pTL4|l : ?^ = 7^(0^^, <) = 7^(0", tt"). 
TC given by Eq. (|III.4p is the function of the gradient vari- 
ables The relation {x,H} = 0, Ea. (|III.10|) . is the 
necessary and sufficient condition to present the Hamil- 
tonian function as a function the projected variables: 

H = H(c,0- 

It is clear that the H is defined finally on the intersec- 
tion of the submanifolds x = and 57 = 0, being thus a 
fimction of the £,s{£,)- Eci- l|III-4|) represents its extension 
to the unconstrained phase space using the full gradient 
projection. 

The statements (i) - (iii) are of the general validity for 
gradient projections. The statement (iii) has been proved 
as such above, the other two ones are proved below. 

D. Dirac bracket calculated by gradient projection 

The phase flow associated to a function g = g(^) de- 
fined on the submanifold Ga = has an ambiguity since 
one can add to g — g{£,) a linear combination of the 
constraints XaGa where Aq are undetermined parameters. 
The phase flow Lg[.] applied to a function / = /(^) suf- 
fers from this ambiguity also: 

iff!/] ={/,5} + Aa{/,^a}. (V.22) 
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The submanifold Qa 



should, however, be invari- 
This equation allows to find 



ant, i.e., LglGa] = 0. 
Aa — Aa(^). Substituting Aa into Ea. ljV.22|l . one gets 
the Dirac bracket 



{f,9}- 

{f,9}D 



{f:Qa}\\{Ga,gb}\\-'{gb,9} 



(V.23) 



where /, g, and Qa are functions of 5. The Dirac bracket 
defines a phase flow generated by a function g = g{£^) 
within the constraint submanifold Qa = 0. 
Using Ea. (|V.ip . one finds at = 



= {/(0(0),5(0(e))}. (V.24) 

This is the analogue of the statement (i) made in the 
previous subsection for the full gradient projection. The 
gradient projection can therefore be used to calculate the 
Dirac bracket. 

There is an analogue of this statement for the partial 
gradient projections also. Let us suppose that the second- 
class constraints Qa = ^ split into the canonical pairs: 
XA = and = such that {xai Xb\ = 0, {^a, = 
0, and det ||{x^,r2B}|| ^ 0. The Dirac bracket becomes 

{/,.9}i3 - {/,.9} (V.25) 

- {!.^a}\\{xa.^b\\\-\xb.9}- 

Let 6(0 and 6(0 be partial gradient projections such 
that Xa(6(0) = and 17a(0(0) = identically (cf. 
Eas. llV.3p and ljV.51) '). If we replace arguments of the 
functions / and g with 6(0 or 6(0j the last two terms 
vanish due to IIV.4|) or (|V.6|I . respectively. The Pois- 
son bracket for the projected variables then coincides 
with the Dirac bracket for the canonical variables ^ con- 
strained to the submanifold xa = or ^a = 0: 

{/(6(0),ff(6(0)}i5 = {/(6(0),ff(6(0)},(v.26) 
{/(6(0),3(6(0)}d = {/(6(0),ff(6(0)}- (V.27) 

The arguments of the functions represent the partial gra- 
dient projections hke in Eas. (|V.9|) - HV.11|) and (jV.lSp - 



Eas. l|V.26(l and l|V.27|l are sufficient to calculate the 
Dirac bracket given that the partial gradient projections 
are constructed. 

This completes the proof of the statement (i) from the 
previous subsection, extended to arbitrary Hamiltonian 
systems. 

Turning to the point (ii), it is sufficient to notice that 
the Poisson bracket, e.g., {6jCs} determines a variation 
of the along the submanifold Qa—^- This submanifold 
is parameterized with the . The Poisson bracket is thus 
a function of the again. The involution relations for 
{6iCs} define therefore an algebra. Similar arguments 
apply for the partial gradient projections. Eas. ljV.9|l - 
(|V.ll|l and ljV.18(l - ljV.20|l represent therefore an illus- 
tration of this statement. 



E. Constraint functions Qa constructed by gradient 
projection 

The statements of Sect. 4 can be proved using the 
gradient projection method. 
The vector fields 



dOa 



(V.28) 



determine phase flows associated to the constraint func- 
tions Qa- These fields are non-singular, i.e., do not 
vanish in a neighborhood of the constraint submanifold 
^ = {C ■ ^a(0 = ya}- The opposite would mean 3a 
such that ljV.28|l vanish at some point ^ G JF. It follows 
then that {Qa,Qb} = V5. This is in contradiction with 



det{Qa,Qb}^0 



(V.29) 



which holds, by assumption, everywhere on and, by 
continuity, in a neighborhood of T. In Ea. (jV.28|) . = 



-lij where 



En 

-En 



(V.30) 



with En being the n x n identity matrix (cf. Ea. (|IV.2|) '). 
In what follows, we denote phase flows ljV.28p briefly as 
IdQa. 

Let Ta = l^- QaiO = 0}. The condition IdQaiO ^ 
e J^a is stronger than the one mentioned above. It 
looks evident, since any constraint hmction ^a(0 in a 
neighborhood of J-a can be redeflned to assign the gradi- 
ent a definite direction. 

is a subspace of the dimension 2n — 1. The inter- 
section of all J-a gives The tangent space to Ta at a 
point denoted as T^Ta, is skew-orthogonal to IdQa{£,). 
Indeed, if dQa = then d^ G T^^Ta, and we obtain 

r^^Ukde = -^dC = 0. (V.31) 

The space T^Ta has the dimension 2?! — 1. Among 
the vectors of T^Ta one can find IdQa{£,), since IdQa{£,) 
is skew-orthogonal to itself. T^Ta is therefore a skew- 
orthogonal complement of IdQa (0 ■ 

One can find b such that {^a(Oj^&(0} 0. As dis- 
cussed above. 



IdQaiO e T^Ta, 

idQbiO e r^^Ff,, 



and furthermore, 



IdQaiO 

IdQbiO 



T^^b, 
T^Ta. 



(V.32) 
(V.33) 

(V.34) 
(V.35) 



From the other side, any vector d^ G T^J-a H T^J-b — 
T^i^a n Tb) is skew-orthogonal to IdQa{£,) and IdQi,{£^). 
As a consequence of Eas. (IV.34(l and l|V.35|l . one gets 



IdQaiO^Wa n Tb), 
IdQbiO i n Tb). 



(V.36) 
(V.37) 
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The subspace Ta^Th and, respectively, T^(Ta'^^h) have 
the dimension 2n — 2. The vectors IdQa{0 ^^nd IdQb{£,) 
are hnearly independent and form a two-dimensional 
space 2?afc(0 which is a skew-orthogonal complement of 
T^i^a n Tb) such that VabiO n T5(J'a n J-fc) = 0. 

Let us consider motion of a particle with a Hamiltonian 
function Qi,- 



In virtue of Eq. ljV.35|l , the phase flow IdQb 
does not lie in the submanifold J-a entirely and therefore 
crosses it. Let ia(C) be time needed for particle located 
at C ^ -^a in a neighborhood of & to cross !Fa at 
some point ?7 ^ ^. The equations of motion look like 



dta 



IdGbiO- 



The derivative of ta(C) along the phase flow IdQb{C) with 
respect to time is, by definition, equal to unity: 



{taicigb{0} = i- 



(V.38) 



One may interpret, equivalently, —GbiC) ^^s a time needed 
to cross the submanifold by a particle located at 
C ^ J-b- The motion of such a particle is described by 
a Hamiltonian function ia(C)- 

The function ta{C) vanishes for ( = rj <E Ta- At any 
point of Ta, dQa and dta vanish for dij e TrfTa- There 
exists only one drj ^ T^jJ-a such that dQa ^ and dta ^ 0. 
It means that Vdr/ dQa is proportional to dta and IdQa is 
in turn proportional to I dta ■ 

The first canonical pair Qa = ta and Qb = Qb is thus 
constructed. 

Let us consider the full gradient projection ^i(^) onto 
the submanifold Ta H J-f,, using the constraint functions 
Qa and Qb. One gets Qa{£.i{0) = 0, Qb{ii{0) = 0, 
whereas the equations Qc{S,i{C)) = for c 7^ a, 6 are sig- 
nificant to determine the location of the constraint sub- 
manifold T, owing to shifts along the phase flows IdQa 
and IdQb- A complete set of equations for T can be taken 
to be 



Qa{0 = 0, 

GbiO = 0, 



for c 7^ a, 6. In virtue of Ea. HV.2p . 

{Qa{O.Qc{Um = 0, 

{Qb{o,Qc{ii{m = 0. 



(V.39) 
(V.40) 
(V.41) 



(V.42) 
(V.43) 



Eas. ljV.40p - ljV.4ip determine T uniquely, so the deter- 
minant of the Poisson bracket relations between the 2m 
constraint functions is not zero. The functions Qa and Qb 
have the vanishing Poisson brackets with the rest ones, 
so 



det{^,(a(C)),a<i(a(C))}^0 



(V.44) 



where c, d take 2m — 2 values (c, d ^ a,b). 

In the remaining set of the constraints, one can flnd 
c,d such that {5c(Ci(C))' and repeat the 



arguments we used earlier. The analogue of Eq. 
looks like 

{ic(C),a<i(a(C))} = l- (V.45) 

The Poisson brackets of the left side of this equation with 
QaiO and Qa{C) vanish. The Jacoby identity yields 

{{tciO,QaiO},Gd{CliO)} - 0, (V.46) 

{{tc{c),Gb{c)},gd{uo)} = 0. (V.47) 

The Poisson brackets of tdC) with Qa{0 and Qb{C) re- 
main therefore constant along the phase flow IdQd{^i{C))- 
At the submanifold = {( : ^c(6(C)) = 0}, IdtdO is 
proportional to IdQc{£,i{C))- Those brackets vanish at J-^, 
and furthermore, vanish for ( ^ T'^. Ea. (jV.2|l suggests 

thent,(C) =tc(ei(0)- 

The second canonical pair QdC) — tc{£.i{0) and 
GdiO = Gd{£,i{0) is thus constructed. 

The proof can be completed by induction. We con- 
sider the full gradient projection ^2(C) onto the subman- 
ifold Ta Ci J-b Ci Tc n Td along the commutative phase 
flows generated by Qa{C), Gb{0, Gc{C), and Qd{C)- These 
constraint functions have the vanishing Poisson brackets 
with the 2m — 4 remaining ones Qe{^2iC)) (e 7^ a,b,c,d). 
The latters constitute a complete non-degenerate set to 
determine the constraint submanifold T uniquely, and so 
on. 

At the end, one gets in a neighborhood of ^ a symplec- 
tic basis (jIV.23|) . 



Hamiltonian Ti. constructed by gradient 
projection 



The Hamiltonian H of Sect, 
with the help of Eg. ljV^ : 



4 can also be constructed 



(V.48) 



where ^s(^) are gradient projections deflned by Qa- 

Ea. ljV.48p and the algorithm described in Sect. 4 give, 
apparently, an equivalent Hamiltonian function, since the 
Hamiltonian flows on the constraint submanifold coin- 
cide. It can be demonstrated by the comparison of the 
Dirac brackets: 



{e,n{mD^{e,niO}D 



(V.49) 



which holds due to Ea. ljV.24p . 

Applications of the gradient projection method to con- 
structing the quantum deformation of the Dirac bracket 
can be found in |40| . 



VI. DIRAC QUANTIZATION OF SPHERICAL 
PENDULUM 

In Sect. 2, we modifled the initial Lagrangian pi.2jl 
on the constraint submanifold x = = to make 
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more transparent the origin of the underlying dilatation 
gauge symmetry. Let us formulate the Hamiltonian dy- 
namics of the spherical pendulum starting directly from 
Lagrangian C + Xhncj). The straightforward follow-up to 
the Dirac's scheme (see also [23,|3^J) leads to the Hamil- 
tonian dynamics described in Sect. 3: 

Using the Dirac's scheme, we obtain the primary 
Hamiltonian Tip — ^tt^ — Alnt/i and the primary con- 
straint ^0 — T^x ~ 0, where tt" are canonical momenta 
associated to the canonical coordinates and 7r> is the 
canonical momentum associated to the Lagrange multi- 
plier A. The canonical Hamiltonian becomes Tic = Tip + 
UTTx where u is an undetermined function of time. The 
secondary constraints Qa+i = {Ga,Ti-c} can be found: 
Gi = \n<j>, 02 = 07r/(^2^ 03 = ^2/(^2 _ 2(07r)V(/.4 + x/cb\ 
The last constraint Qs = allows to fix A in terms of 0" 
and tt", no new constraints then appear. 

The dimension of the phase space can be reduced by 
eliminating the canonically conjugate pair (AjTTa). We 
solve equations Go = with respect to ttx and Ga — 
with respect to A and substitute solutions into Tic- The 
result is 

H'^ = ^Tr^ + {Tr'^ -2{(j)Tr)y(j>^)ln(P. (VI.l) 

There remain two constraint functions Gi and G2, such 
that {Gi,G2} = 1/0^- These are the second-class con- 
straints. The Hamiltonian H'^ is first-class: {5i,7i^} = 
G2 - Gi{Gi,M and {G2,'H'^} = -^i{^2,A} where A is 
determined hy G3 — 0- 

In order to split the constraints into the gauge-fixing 
condition and the gauge generator, we should construct, 
according to the previous Section, Ga and Ti. starting from 
Ga and H'^. 

In Sect. 3, we already had the same constraint sub- 
manifold described by functions x = In ^ and Q = (/jtt 
satisfying {Xj f^} — li so one can set Gi = ln0 and 
G2 = ^Ti". 

The Hamiltonian Ti. can be constructed with the help of 
Eq.lfOsl) . Let us combine Eqs.Ejj), fV^S)! and (fyTTG)! . 
(|V.17|I to get the full gradient projection: 

c = r/ci^, (VI.2) 

TT^ = (/)7r" - (/)"07r/(/). (VI.3) 

Replacing 0" and tt" in by variables (j)" and irf, re- 
spectively, we get 

n = \t:1 (VI.4) 

which coincides with the Hamiltonian (|111.4|l . 

The difference Ti. — Ti'^ is of the second order in the 
constraint functions. This is sufficient to have identical 
Hamiltonian flows on the constraint submanifold. 

The equivalence is thus demonstrated, so further dis- 
cussion can go in the parallel with Sect. 3. In the exam- 
ple considered, the constraint functions are of the second 
class, whereas in Sect. 3 they appear as the gauge-fixing 



condition, x = and the gauge generator, il. This sug- 
gests that the interpretation of second-class constraints 
of a holonomic system is a matter of convention. 

VII. SECOND-CLASS CONSTRAINTS AS 
GENERATORS OF GAUGE SYMMETRIES AND 
GAUGE-FIXING CONDITIONS 

The type of constraints appearing within the Hamilto- 
nian framework depends on the form of the corresponding 
Lagrangian. In case of the spherical pendulum, starting 
from Ea. (|II.2p one arrives to the constraints x = and 

= as to the second-class constraints of the Hamilto- 
nian framework. Starting from Ea. ljII.ll|l . the constraint 
X = appears as a gauge-fixing condition, whereas the 
constraint = as a gauge generator of the symme- 
try group. Lagrangians j!L2)l and pLTT|l are equivalent 
at the classical level, so they lead to the same classical 
dynamics. 

It is possible, therefore, at least in the case of spheri- 
cal pendulum, to interpret second-class constraints as a 
gauge-fixing condition and a gauge generator, and vice 
versa. 

We wish to discuss whether this is a common situation. 
Any set of admissible gauge-fixing conditions and gauge 
generators can be treated as second-class constraints. 
This statement is widely used for quantization of gauge 
theories (see, e.g., [lT|l. 

The reverse statement represents apparent interest 
also. The second-class constraints systems are anal yzed 
from this point of view by Mitra and Rajaraman |29l |. 
In this Section, we discuss additional details about the 
underlying gauge symmetries of second-class constraints 
systems. 

A. Converting a gauge system into a second-class 
constraints system 

Let a primary Hamiltonian ?i of a system be gauge 
invariant. The generators of gauge transformations fi^i = 
are such that 

{nA,H} = Qb, (viLi) 
{nA,nB} = nc. (vil2) 

The gauge-fixing conditions xa = fulfill 

{xa.Xb} = 0, (VIL3) 
det{xA,nB} ^ 0. (VIL4) 

After imposing gauge-fixing conditions, the system be- 
comes equivalent to a system described by a first-class 
Hamiltonian Ti' defined on the submanifold of second- 
class constraints Ga — {xa,^a)- The Hamiltonian Ti' 
can be constructed from the canonical Hamiltonian of 
the original system. Tic = Ti + Xb^b, requiring 

{XA, He} = {XA, n} + Xb{xa, flBj^O. (VIL5) 
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The gauge parameters As can be fixed from this equation 
in terms of the canonical variables. We thus get the first- 
class Hamiltonian H.' = He and second-class constraints 
Ga satisfying {^/a, W} « and det{ga,Gb} 7^ 0- 

The gauge-fixing is equivalent to converting the orig- 
inal system into an equivalent one described by a first- 
class Hamiltonian and second-class constraints. 



B. Converting a second-class constraints system 
into a gauge system 

It was shown by Dirac 0,01 that first-class constraints 
imply the presence of unphysical degrees of freedom the 
evolution of which is not fixed by the Hamilton's equa- 
tions. The dynamics is self-consistent provided the un- 
fixed degrees of freedom do not belong to the set of phys- 
ical observables. This is the case of gauge-invariant sys- 
tems which have gauge degrees of freedom. According to 
the Dirac's constraint dynamics, having first-class con- 
straints is a necessary condition for the existence of a 
gauge symmetry. The reverse statement is less obvious: 

Gauge transformations of the canonical coordinates 
have the form 

6r = {r,^A}0A (V11.6) 

where 9a are infinitesimal functions of time. Since ^Ia 
depend on the canonical coordinates and momenta, the 
variations 5(j)°' depend on the coordinates and momenta 
also, unless fJ^ are first degree polynomials of the canon- 
ical momenta. Subject to a Legendre transform, 50" be- 
come, in general, functions of the coordinates and veloci- 
ties. The gauge theories such as QED and QCD have the 
constraint functions Qa as the first degree polynomials 
of the canonical momenta. In such cases Scj)" = 6(j)"'{4>^) 
do not depend on velocities and define thereby trans- 
formations on the configuration space A4 = {4>°') which 
induce, in turn, transformations in the tangent bundle 
TM = (0", 

The first-class constraints systems correspond to more 
general class of gauge-invariant systems with 50" — 
S(f>°'{(j)^ ,(f>^) and, probably, new auxiliary variables. In 
what follows, we distinguish thereby between gauge 
transformations of the form dcj)" — 54)"'{(t)^) and gen- 
eralized gauge transformations with velocity dependent 
parameters and/or new auxiliary variables. 

Having a first-class constraint system is a necessary, 
but not a sufficient condition for an equivalent gauge- 
invariant system to exist in the original configuration 
space. 

The possibility of constructing gauge- invariant systems 
in the unconstrained phase space, equivalent to second- 
class constraints systems upon a gauge- fixing, has been 
analyzed by Mitra and Rajaram an |29|| . 

The applications discussed [Sfl |^ have constraint 
functions, associated to gauge transformations, as poly- 
nomials of the canonical momenta of the degree less or 



equal to unity. The gauge transformations are not veloc- 
ity dependent, although involve new auxiliary variables. 
The gauged systems require generally an extended con- 
figuration space. 

// global symplectic basis in the space of constraint 
functions can be found, equivalent generalized gauge sys- 
tems for second-class constraints systems can be con- 
structed in the unconstrained phase space. 

By passing over to the standard canonical frame, one 
can always select ^a as first canonical momenta. This is, 
however, not sufficient to have gauge invariance on the 
physical configuration space. Before doing a Legendre 
transform, one has to pass first to a canonical coordinate 
system where the coordinates 0" constitute a physical 
configuration space. This will mix up the canonical co- 
ordinates and momenta, preventing from having ^Ia as 
the first degree polynomials. 

It is worthwhile to notice that gauge invariant observ- 
ables are not measurable if they involve auxiliary de- 
grees of freedom. The sums of the vector potentials of 
the massive electrodynamics and the derivatives of the 
Stiickelberg scalar are gauge invariant. They do not be- 
long, however, to the set of physical observables. The 
equivalence with the ordinary gauge systems, where sets 
of gauge invariant quantities and physical observables co- 
incide, is therefore not complete. 



C. Gauge-invariant systems as holonomic systems 

The quantization of gauge theories which appear under 
first-class constraints is studied in many details (see, e.g., 
Refs. 0, 0). The second-class constraints systems 
are usually quantized by employing the BFT formalism 
[Tol IT^ Eof that converts second-class constraints into 
a first class by extending the original phase space to a 
higher dimension. 

The BFT algorithm combined with the Fradkin- 
Vilkovisky quantization scheme Jj if applied to a sys- 
tem in a 2n-dimensional phase space ends up with an 
extended phase space of a dimension larger or equal to 
2n-\-12m, since the 2m constraints convert into first-class 
constraints, each of which requires an independent gauge- 
fixing condition. The remaining at least 8m degrees of 
freedom appear as ghost variables. 

The problem of constructing gauge-invariant systems 
equivalent to second-class constraints systems in the orig- 
inal phase space is discussed in Ref. 29] . The existence 
of the global symplectic basis for the constraint func- 
tions is specified as sufficient conditions for the existence 
of gauged partners associated to the second-class con- 
straints systems. The method '291, being successful in 
removing auxiliary fields from the original phase space, 
generates auxiliary fields in the effective Lagrangians. 

The systems of pointlike particles under holonomic 
constraints have the natural gauge counterparts both in 
the phase space and configuration space |27| . 

In holonomic systems, second-class constraints split 



into gauge-fixing conditions, xa = 0, and gauge genera- 
tors, flA- Such systems can be quantized further as gauge 
theories. The BFT method, if apphed, would result to an 
extended phase space of the dimension 2n + 6m, contain- 
ing 2m Lagrange multipliers and 4m ghost variables. The 
underlying gauge symmetry allows to reduce the number 
of auxiliary fields within the quantization scheme 0, Q . 

D. Frobenius' condition for holonomic constraints 

The holonomic constraints are defined by a constraint 
submanifold TV C in the configuration space Ai. It 
comes then out automatically that particle velocities be- 
long to the tangent space TA/" of the constraint subman- 
ifold. 

An n — m-dimensional tangent subspace TJ\f C TA4 
of an n-dimensional configuration space A4 can be de- 
fined in general through m covector fields UAa with 
A = l,...,m and a — l,...,n by imposing constraints 

UJAar = (VII.7) 

which are velocity dependent and therefore non- 
holonomic. 

The Frobenius' integrability conditions (see, e.g., [28|'). 




(0?<^^ - <^^</)f ) = 0, (VII.8) 



for tangent vectors 0" and 02 satisfying l|VII.7|l . if ful- 
filled, implies the existence of a submanifold TV the tan- 
gent space of which coincides with the tangent subspace 
(|V11.7|I . In such a case, the constraints (|VII.7() can be 
replaced with holonomic constraints XA = 0, identifying 
the constraint submanifold J\f d M, without modifying 
the dynamics. 

Eas. (|VII.8|l specify non- holonomic systems which have 
gauged counterparts in the original configuration space. 

For an n = 3 Euclidean space, Eas. (|VII.7p define, 
for m = 1, a submanifold orthogonal to the vector uja- 
Eas. HVII.8p teU us that the vector field cJ is a potential 
field provided that rotuj = 0. One can find a potential 
function x such that to — Vx- The sets of x = constant 
define, in our case, possible constraint submanifolds of 
an equivalent holonomic system. 

The tangent subspace determined by Ea. (|II.3p coin- 
cides with the tangent space of S^~^. The covector field 
ijJa = satisfies the Frobenius' condition. 

VIII. SUPPLEMENTARY CONDITIONS FOR 
WIGNER FUNCTIONS 

The splitting of second-class constraints into con- 
straints associated to gauge transformations and gauge- 
fixing conditions is not unique. Every pair {xa,^a) can 
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be transformed, e.g., as xa ^a, ^a ~Xa- The 
symplectic group Sp(2m) of transformations mixes xa 
and VLa further, keeping the Poisson brackets invariant. 
Canonical transformations, furthermore, do not modify 
the Poisson bracket of the constraints also. There exists 
therefore an affluent spectrum of representations and we 
have to demonstrate that the physical content of theory 
is invariant with respect to the choice of representation. 

In the general case of second-class constraints, we do 
not have any criterion how to discriminate between first- 
class constraints associated to gauge transformations and 
gauge-fixing conditions. We show, conversely, that the 
way the second-class constraints are split is not impor- 
tant for quantization. 

A useful hint towards that conclusion is delivered by 
the path integral representation for the evolution opera- 
tor ljl.6|l which is invariant with respect to linear trans- 
formations of Qa- 

The problem is to demonstrate that the supplementary 
conditions VIa"^ — can be replaced by xa^ = 0. In case 
of the spherical pendulum, this is almost evident, since 
Cl'^ = means 4* ^ / (p) whereas = means 

^^f ^ 5(x)*((/)") = (5(x)^'(0"/(/>). The essential depen- 
dence comes from the angular part of the wave function, 
while the radial part S{x) is factorized, being commuta- 
tive, Ea. llIII.10|l . with the Hamiltonian and the S'-matrix. 
It can be absorbed thereby by an overall normalization 
factor. 

These arguments can apparently be extended to an 
arbitrary case. According to the Dirac's prescription, the 
physical wave functions are annihilated by the constraint 
operators associated to gauge transformations 

17^* = 0. (VIII.l) 

To simplify the matter, we pass to the standard canon- 
ical coordinate system and choose the constraint func- 
tions XA as first m canonical coordinates qa and the 
constraint functions Qa as first m canonical momenta 
PA- The remaining canonical variables are {q*,p*) con- 
stitute the physical phase space p^("~'")_ Ea. HVIII.ip 
gives then, in the coordinate representation, a wave func- 
tion of the form ^ — '^{q*). Such a wave function 
has an infinite norm, since the integral / (i"g|^'(g*)|^ 
diverges. In the momentum representation, we would 
get *(p) = (27rn)™nA'5(P^)*(P*)- has an infi- 

nite norm (27r?i)™5™(0) provided ^'(p*) is normalized to 
unity. An infinite but numerical factor can be included 
into the norm of ^'(p). 

Let us check whether wave functions satisfying 

XA* = 0. (VIIL2) 

have a physical sense. In the coordinate representa- 
tion, one gets * — Y\A^ilA)'^{'l*) and in the momen- 
tum space ^' = vE'(p*). These states have infinite norms. 
There is an apparent symmetry p q, q ^ —p between 
wave functions satisfying ljVIII.l|) and (|VIII.2|I . 
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Owing to the factor YIa ^(qa), conditions (|VIII.1|) and 
(|Vlii.2() single out the same set of functions, the nontriv- 
ial dependence of which comes from the physical variables 
q* (p*) only. 

The Wigner functions |37l | provide complete informa- 
tion on quantum systems. The association rules of op- 
erators in the Hilbert space and functions in the phase 
space are discussed a long time 36^ 37, 38, 3^ The 
quantum mechanics can be reformulated using the Groe- 
newold star-product [s^ representing a deformation of 
the usual pointwise product of functions in the phase 
space. The Wigner functions and constraint functions are 
defined in the phase space, so it is natural to discuss sup- 
plementary conditions in terms of the Wigner functions. 
In this approach, the equivalence between Eas. ljVIII.l|) 
and ljVIII.2p becomes more transparent. In addition, 
more general group of supplementary conditions, equiv- 
alent to Eas. HVIII.il) and (|VIII.2|) . can be formulated in 
terms of the Wigner functions. 



A. Probability density localized on the constraint 
submanifold 

Let us start from the standard canonical frame [s^l 
where the constraint functions XA are the first m canon- 
ical coordinates qa and the constraint functions Qa are 
the first m canonical momenta pA- The probability den- 
sity in the phase space F^" can be written as follows: 

W{q,p) = {2nhy^l[S{qA)S{pAW*iq*,P*)- (VIII.3) 
A 

Identifying the W{q,p) with the Wigner function in the 
unconstrained phase space and using the Wigner trans- 
form, one gets the density matrix 



and p replaced with UpU+ , Eas. (|VIII.5p and (|VIII.6p be- 
come 



Yl6C^^^)p.{q\q*'). (VIIL4) 



GaP + pQa = 0. 



(VIII.8 



Eas. (|VIII.8|l are necessary and sufficient conditions to 
have the representation (|VIII.4(I in the standard canoni- 
cal frame. 

The supplementary conditions f VIII . ig)) cannot he for- 
mulated in terms of a wave function, since the density 
matrix in the unconstrained configuration space does not 
correspond to a pure state, even if system on the con- 
straint submanifold is in a pure state. 

In an arbitrary frame and in the classical limit, 
Ea. (|VIII.3p looks like 

w[0 = {2^hY''\{5{ga)Vd^t{gZQb}w.[UO)- 

(VIII.9) 

Note that Ha ^i^a) acts as a projection operator, so that 

Y{a5{Qa)f{0 ^Y{aKQa)f{UO) ^ f iO ■ 

In the classical limit, the Wigner function satisfies 



(VIII. 10) 



The complete phase-space analogue of quantum 
Eas. ^vnm can be formulated in terms of the symmetric 
part of the Grocncwold star-product Isl] to give 

SaiO ° W{0 = 0. (VIII.ll) 

The star-product has the following decomposition: 

/(O * 9(0 = /(O ° 5(0 + y /(O A cM) (VIII.12) 



where 



/(O 0.9(0 = /(0cos(-P)5(e), (VIII.13) 
/(OA 5(C) - /(0|sin(j7')g(0 (VIII.14) 



and 



It satisfies the operator equations 

qAp + pqA = 0, 
PAP + PPA = 0. 



(VIII.5) 
(VIIL6) 



The first equation implies p S{qA + i'a)- The second 
one implies that p does not depend on qA — q'a- 

The density matrix p*{q*,q*') is normalized to unity, 

so 



/ 



p{q,q)d"q=l. 



(VIII.7) 



One can make a unitary transformation to pass to an 
arbitrary set of operators associated to the canonical vari- 
ables. In terms of 

Ga = U{qA,PA)U^ , 



V 



-P- 



d d 



is the so-called Poisson operator. In the limit h 



lim/(OA.g(0 



{/(0,5(0}- 



0, 



(VIII.15) 



The Poisson bracket {/(Oi 5(0} coincides generally with 
a function associated to the commutator —i/h[f,g] to 
the lowest order in the Planck's constant only. The 
skew-symmetric part of the Groenewold product provides 
a generalization of the Poisson bracket which is skew- 
symmetric with respect to two functions, real for real 
functions, coincides with the Poisson bracket to the low- 
est order in the Planck's constant, satisfies the Jacoby 
identity, and keeps the association rule for commutators. 
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This part of the Groenewold product is known also as 
the "sine bracket" or Moyal bracket 39] . 

In the Umit h ^ 0, the quantum condition l|VIII.ll(l 
recovers the classical one (|VIIL10|I . These two conditions 
coincide in the standard canonical frame. 

The normalization condition Ea. (|VIII.7p holds in the 
quantum case. The Wigner function satisfies 



w{0 



{27rhy' 



1. 



(VIII.IC 



The classical limit of the Wiener function Ea. HVIII.9p 
provides the conventional normalization on the con- 
straint submanifold for VF,(^s(^)). 

The Groenewold star-product is non-local. The usual 
pointwise product of two functions like in Ea. ljVIII.10|) 
has no quantum counterpart. It can be treated as a clas- 
sical limit of a quantum operator equation only. By con- 
trary, Ea. (|VIII.lip is the phase space analogue of the 
quantum operator equation (|VIII.8|) . 

Eq. I|VIII.3|) accomplishes a trivial extrapolation of the 
Wigner function from the constraint submanifold: The 
density is set equal to zero when ^ does not belong to 
the constraint submanifold. The Wigner function is not 
a smooth function across the constraint submanifold, so 
it is hard to formulate using this approach an evolution 
equation similar to the Liouville equation in the uncon- 
strained phase space. 



B. Probability density localized on and outside of 
the constraint submanifold 



One can require that at any given point ^ the density 
be the same as at = Cs(C)- The gradient projections S,s 
can be constructed, as discussed in Sect. 5, using phase 
flows generated by the constraint functions Qa to solve 
equations Ga{£,s{0) 0^ 

The analogue of Eq. (|VIII.3ll reads 



(VIII. 17) 



Such a density has infinite norm, since there are direc- 
tions in the phase space F^", crossing the submanifold 
Qa = 0, along which the density remains constant. These 
directions are determined by the constraint functions. 
Ea.fOl tells 



(VIII.18) 



In the standard canonical frame, where qa — Xa and 
PA = ^A, the gradient projections can be easily con- 
structed: 



6(0 = (o,-z*,o,p*) 



(VIII. 19) 



where ^ — {qA,q* ,Pa,P*)- Ea. llVIII.17p then simplifies 
to 



Note that 

= (0,g*,PA,P*), (VIII.21) 
= iqA,q*,0,p*), (VIII.22) 

If we apply the Wigner transform, we get the density 
matrix 

= l[S{qA-q'A)p*{q\q*')- (VIII.23) 



It satisfies 



or, equivalently. 



qAP - pqA = 0, 
PAP ~ PPA = 



GaP - pGa = 0. 



(VIII.24) 
(VIII. 25) 



(VIII. 26) 



The phase space analogue of these operator equations 
looks like 



GaiO^W{0=0. 



(VIII.27) 



This condition is in agreement with its classical counter- 
part E q.(|VIII.18|) by virtue of Eg. lVIII.lSp. 

Eas. (|VIII.26|) are distinct from Eas. HVIII.8|) . This 
is a consequence of different extrapolation schemes 
of the density to the unconstrained phase space. 
Eas. ljVIII.26p arc the necessary and sufficient conditions 
for Ea. jViil.23t . 

The dependence of the density matrix on (qA, 9^) does 
not factorize, although — p for pi = p*. The latter, 
even when fulfilled, is not sufficient to have a pure state. 
Ea. (|VIII.23p describes, in particular, a noncoherent sum 
of pure states with different momenta pA- The system 
is thereby identified as a mixed state. It cannot be de- 
scribed by a wave function. The constraint equations 
(jVIIL26p . respectively, cannot be formulated in terms of 
a wave function '^{q). 

If one works with density matrices in the represen- 
tation ljVIII.3p or ljVIII.20p . it is not important how 
the second-class constraints were split. Eas. (|VIII.8p and 
(|VllI.26p are symmetric explicitly with respect to the 
interchange XA ^A, ^A ~XA, linear transforma- 
tions of the constraint functions, and furthermore, with 
respect to unitary transformations in the Hilbert space. 
In the classical limit, they are invariant with respect to 
canonical transformations. 

The normalization condition of the Wigner function 
involves a projection operator 



W{q,p) = W.,{q*,p*). 



(VIII.20) 



-i2m 



. 2m 



(VIII.28) 
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in terms of which the norm is calculated as follows l4C 



(VIII.29) 



where is the Weyl's symbol of the operator 

(IVIIL28II . Note that *p is commutative with p in virtue 
of ^VTTT.2611 . 

The Wigner function appears as a smooth function, 
so it is an appropriate object to describe an evolution 
of the system on line with the Liouville equation in the 
unconstrained phase space. A quantum extension of the 
the Liouville equation for W{(,) satisfying Eas. (|VIII.27|) 
is given in Rcf . ,4QJ . 



C. Probability density with mixed localization 

Let us keep, e.g., canonical momenta on the constraint 
submanifold and use the gradient projection to extend 
the density as a function of the canonical coordinates 
away from the constraint submanifold. The Wigner func- 
tion is given then by 



W{q,p)^{2TThr\{5{pA)WM\ 



P 



(VIIL30) 



Applying the Wigner transform, we get the density ma- 
trix 



which satisfies 



PAP 
PPA 



0, 




or, equivalently, 



QaP = 0, 
pClA = 0. 



(VIIL31) 



(Vin.32) 

(VIIL33) 



(VIIL34) 
(VIIL35) 



These equations provide the necessary and sufficient 
conditions for Ea. (|VIII.3ip . No dependence on the 
{qA,q'A) appears. If pi = in the coordinate space q* on 
the constraint submanifold then = p in the coordinate 
space q (owing to an infinite normalization factor). 

The system is allowed to appear in a pure state and 
can be described by a wave function VE'(q) accordingly. 
Eqs. \VIII.34\) and f VIII. 3^) can be reformulated in terms 
of wave functions to match the Dirac's supplementary 
condition Eo. iVIII.l\) . 

The mixed localization scheme breaks the symmetry 
XA — * ^A, ^A —XA from the outset. However, it al- 
lows to work with wave functions. The other methods we 
discussed lead, in the unconstrained configuration space, 
to mixed states. 



The quantum analogue of Eas. (jVIII.34|) and HVin.35|l 

in the phase space can be formulated in terms of the 
star-product: 



nAiO^wiO = 0, 
w{0*i^A{0 = 0. 



(VIIL36) 
(VIIL37) 



As a consequence of Eas. HVIIL34|) and HVIIL35|), we 
have {f2A,W^} = in the classical limit. Ea. HV.4|l 
tells then that W{^) = p{£,u{£,))- The second classi- 
cal condition il^W^ = results in the Wigner func- 
tion W{^) = Yl S{flAiO)P*{^v{0) ■ Combining these two 
equations, we obtain 

W{0 - {27rhrl[S{^A{0)W4U0), (VIIL38) 



which is in agreement with Ea. (|VIII.30p . We used here 
V.a{^u{0) = ^a{0 which is valid up to the second order 
in the constraint functions XA- 

It is possible to establish a relationship with the results 
of the previous subsection. Indeed, one can check that 



VnPs, 



(VIII.39) 



with ps being the density matrix from the previous sub- 
section, satisfies Eas. (|VIII.34|) and (|V11I.35|I . The oper- 
ator ?Pf2 is defined by 



%2 



d"'X 



J]^ exp(^riAAA) 

A=l 



(VIII.40) 



In order to get the normalization condition for W^(^) 
one has to factorize p according to (|VIII.39p . construct 
the Weyl's symbol of ps, and apply Ea. (|VIII.29|l . 

Let us consider the opposite localization: 



W{q,p) ^\[5{qA)W,{q\p*). 



(VIII.41) 



The first m canonical coordinates are kept on the con- 
straint submanifold, whereas the first m canonical mo- 
menta are projected. 

The density matrix has the form 

p{q.q')^J\5{qA)6{q'A)p*{q\q*') (VIII.42) 



It satisfies qAp = pqA = or, equivalently. 



XAP 
PXA 



0, 
0. 



(VIII.43) 
(VIII.44) 



These conditions are the necessary and sufficient condi- 
tions to have the density matrix of the form HVIII. 42|l . 

The system is allowed to appear in a pure state and 
can be described by a wave function. The constraints 
(|VIII.43ll and (|VIII.44p can be reformulated in terms of 
the wave functions to give Ea. (jVIII.2p . 
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As a consequence of Eas. (|VIII.43|l and HVIII.44|) . one 
gets in the classical limit 

WiO =l[^iXAmW4U0l (VIII.45) 



in agreement with Ea. (|VIII.41(l . We used here the rela- 
tion Xa{S,v{0) = Xa{0 which is valid up to the second 
order in f2^. 

The quantum analogue of Eas. (|VIII.43|) and HVIII.44|) 

is given by 

Xa{0*W{0 = 0, (VIII.46) 
W{0*Xa{0 - 0. (VIII.47) 

There exists a relationship with the Wigner function 
of the previous subsection. One can check that 



(VIII. 



satisfies Eas. HVIII.43ll and (|VIII.44|I . The operator is 
defined by 



XaV^ 



(VIII.49) 



in the 



Note that ^ = ^nVx and ilA^n 
operator sense. 

The normalization of the Wigner function in arbitrary 
canonical coordinate system is quite involved: 

Given a wave function in the unconstrained configura- 
tion space, one should construct the density matrix, fac- 
torize it according to Egs. f VIII. 39\) or \VIII.4^ and in- 
tegrate the Wigner function associated to ps in order to 
extract the norm according to Eg. f VII1.29\) . 

Respectively, the Hermitian product of wave functions 
is calculated using the off-diagonal Wigner function and 
integrating it over the unconstrained phase space accord- 
ing to the same prescription. 



original phase space as pure states (mixed localization). 
One of them can be constructed using the Dirac's pre- 
scription l|VIII.ip . The dual condition ljVIII.2p was found 
to be possible also. If we work with density matrices or 
the Wigner functions, one arrives at conditions (|VIII.8|) 
or HVIII. 18() . both are symmetric under the permuta- 
tions xa ^A, ^A — > ~Xa- The constraints imposed 
on the physical states do not depend on the splitting 
Ga = (xAi^a)- 



IX. THE 0{n) NON-LINEAR SIGMA MODEL 

The 0{n) non-linear sigma model represents the field 
theory analogue of the spherical n — 1-dimensional pen- 
dulum. The n = 4 case corresponds to the chiral non- 
linear sigma model due to the isomorphism of algebras 
su{2)l © su{2)r ~ so(4). 

In Sect. 2, we started from the tangent bundle TA4 = 
((/)", 0") of the dimension 2n, defined over the configura- 
tion space A4 ~ ('/'")■ Lagrangian (|II.2|I depends on n 
velocities 0" . On the constraint submanifold l|II.l|) , it de- 
pends on ri — 1 velocities A°'^(t)^ tangent to the constraint 
submanifold. The dynamics of turns out to be inde- 
pendent on other variables. The constraint (jll.l|l reduces 
the effective number of the degrees of freedom to n — 1. 
Thus, a 2n — 2 dimensional tangent bundle over the n — 1 
configuration space can be constructed. It is described, 
e.g., by coordinates analogous to the angular coordi- 
nates in three dimensional space. Lagrangian (|11.11() is 
not degenerate with respect to The coordinates 
constitute the physical configuration space. 

The path integral for the evolution operator in terms 
of the angular coordinates can be treated as a refer- 
ence point for comparison to more involved quantization 
methods. 



D. Discussion 



A. Path integral for the 0{n) non-linear sigma 
model 



In the classical constraint systems, physical quanti- 
ties depend on the probability densities localized on the 
constraint submanifold only. In the standard canon- 
ical frame, the Wigner functions of quantum systems 
are localized on the constraint submanifolds also. This 
is mandatory, since any unconstrained system can be 
treated as a constrained system with constraints imposed 
to remove the added unphysical degrees of freedom. In 
doing so, the unphysical degrees of freedom should not 
modify dynamics of the initial system. This is achieved 
by attributing physical sense to the Wigner functions on 
the constraint submanifold. How to extrapolate Wigner 
functions from the constraint submanifold into the un- 
constrained phase space is a matter of convention. 

We discussed the most evident extrapolations. Among 
them are those which allow to describe systems in the 



Let us construct the path integral for the 0{n) non- 
linear sigma model in terms of the angular variables 
i.e., by solving the constraint equations from the outset 
and compare it with the expression of Sect. 2 derived 
using the underlying gauge symmetry of the 0{n) non- 
linear sigma model. 

The field theory extension of the spherical pendulum 
problem is rather straightforward. In what follows, the 
kinematic variables are functions of = (i,x). 



Let us construct a basis 



a a 



d 

w 

— 9ij J 



(IX.l) 

(IX.2) 
(IX.3) 
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where gij = gij{'&) is an induced metric tensor on the 
submanifold 4'{x) = 1, detl|(7y|| ^ 0. In terms of the 
coordinates -d^ , the field theory extension of Lagrangian 
(|II.11|I takes the form 



(IX.4) 



The Legendre transformation of is well defined. In 
term of the canonical momenta 



9£* 



the Hamiltonian density can be found to be 



(IX.5) 



(IX.6) 



where a — 1,2,3. The non- vanishing Poisson bracket for 
the canonical coordinates and momenta has the form 



B. Pion field parameterization in chiral sigma 
model 



The n = 4 case is especially interesting since it corre- 
sponds to the chiral non-linear sigma model. For n = 4, 
the angular coordinates are defined by 

0" = {cos ^p, sin ^px {cos 9, sin 9 X {cos If, sin tp))) (IX. 12) 

where -d^ = ip, i^"^ = 0, and -d^ = (p. 

The angular distance, 9, between two vectors 0" and 
(/)'" is defined by scalar product cos0 — (p(j>' . The dis- 
tance element becomes 



de^ = dip"^ + sin^ il;{d9^ + i 



'■)■ 



(IX. 13) 



The components of the metric gij can be found using the 
expansion dQ^ = gijdd^dd^ or directly from Eas. (|IX.l|) 
and ljIX.2|) . The Lagrange measure of the path integral 
becomes 



K(i,x),e,(t,x')}-<5;-<5(x-x'). 



(IX.7) 



The Poisson bracket relations for coordinates and 
momenta tt" associated to the tangent velocities vr" — 
Qig^''Sj agree with those discussed in Sect. 4. 

The path integral in the phase space (i9% ft) is given 

by 



Z = 



(IX.. 



The Liouville measure dd^dgi is consistent with 
Ea. (|IX.7l) . The integral over the canonical momenta in 
Eq. (|IX.8I) has a Gaussian form and can be calculated ex- 
plicitly: 



Z = 



nVdet 



\9ij 



|d"-iz?exp<j ^ / dx^C, 



(IX.9) 

The value i/det ^i? gives volume of the config- 

uration space, defined by the metric tensor gij. This 
measure is invariant under the 0{n) group. 

The ^-matrix ljIX.9p can be written in an explicitly 
covariant form with respect to the 0{n) rotations and the 
Lorentz transformations. First, we rewrite Lagrangian 
density (|IX.4|) in terms of the coordinates </>" 



1 



(IX.IO) 



and, second, rewrite the Lagrange measure 



^detU,\\d-~',9 = ^{dx/drfS{x)d-(P (IX.ll) 

where x — In (p- The right side is the same for all func- 
tions vanishing at </> = 1. 

In this form we recover the result of Sect. 2. 



^Jdet\\gij\\d^'d = sin^ ip sin 9 d%pd9d(p 
sin? -0 



-dV, 



(IX. 14) 



with dV = tjj^ dip sin 9 dO dip being an element of the Eu- 
clidean volume. The corresponding Lagrangian can be 
found from Ea. ljIX.4p to give 

1 



M^PY 



sin^ ip 



{{d^9Y + sin' 9{d^^Y). (IX.15) 



The quantization of the chiral sigma model is made 
using an oscillator basis by expanding the around i?* = 
breaking thereby the 0(4) symmetry down to its 0(3) 
subgroup. It can be successful provided that measure of 
the coordinate space is such that det \ \gij\\ = 1. The path 
integrals convert then to the Gaussian integrals which 
can be calculated. The parameterization preserving the 
0(3) symmetry and satisfying the above requirement is, 
apparently, unique. One should rescale the "radius" ip 
according to 



sin^ ipdip = uj^duj. 
This elementary equation gives 



-(-0 — sin ■0 COS ■0) 



1/3 



(IX. 16) 



(IX.17) 



Lagrangian C■^, then becomes 

4 • 2 / 

sm ip 



-{d^ojY + 



{{d^9Y+sin^9{d^^Y) 



2sin'tP'~''"' ' 2 

(IX.18) 

where '0 is a function of ui Ea. l)IX.17(l . The mass term 
breaking the 0(4) symmetry looks like 



Cm = M^{cos'iJj - 1). 



(IX. 19) 
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The quadratic part of the Lagrangian used for the per- 
turbation expansion can be selected as follows 

In terms of the pion fields 

tt" = ti> X (cos 6*, sin 6* X (cos (p, sin (/s)), (IX. 20) 
it takes the standard form 

(IX.21) 

whereas the Lagrange measure is simply the Euclidean 
volume (P-K. The difference (5£* = £* + Cm — -Cp' can 
be considered as a perturbation. 

The pion fields can be parameterized in various ways. 
The problem of ambiguities of the transition amplitudes, 
connected to the arbitrariness of that choice, was dis- 
cussed first in Refs.[4ll |43|. It was shown that on-shell 
amplitudes do not depend on the choice of physical vari- 
ables. This statement is known as the "equivalence the- 
orem" . The method proposed by Gasser and Leutwyller 
1^3 associates the QCD Green functions to amplitudes of 
the effective chiral Lagrangian. Using this method, the 
QCD on- and off-shell amplitudes can be calculated in a 
way independent on the parameterization. 

The S'-matrix is invariant with respect to a symmetry 
group, if both the action functional and the Lagrange 
measure entering the path integral over canonical coordi- 
nates are invariant. The Liouville measure entering the 
path integral over canonical coordinates and momenta 
is always "flat", since we work in canonical basis. The 
quantization gives a non-trivial Lagrange measure, how- 
ever (cf. Ea. ljIX.9|) '). In case of the 0(n) non- linear sigma 
model, there is only one parameterization which makes 
the Lagrange measure flat, i.e., det||f;ij|| — 1. This re- 
quirement is useful for development of perturbation the- 
ory which uses an oscillator basis to convert path inte- 
grals into the Gaussian form. 

The weight factor can always be exponentiated to gen- 
erate an effective Lagrangian 5Ch, in which case x" — <f>°' 
provides desired parameterization also. The exponential 
parameterization of the pion matrix 

[/((^") = 6^^°"^° (IX.22) 

gives, in particular, 

5Ch = -^Msin^UU-'^) (IX.23) 

where a is a lattice size. 5Ch diverges in the continuum 
limit. The non- linear sigma model is not a renormalizable 
theory, so divergences cannot be absorbed into redefini- 
tion of F and M^. Using the mean filed (MF) approxima- 
tion, it is usually possible to keep renormalizations finite. 
The exponentiation of a variable weight factor breaks, in 



general, selfconsistency of the MF approximation of the 
non-linear sigma model. 

Divergences arising from 5Ch could, however, be 
compensated by divergences coming from higher orders 
ChPT loops. From this point of view it looks natu- 
rally to attribute 5Ch to higher orders ChPT loop ex- 
pansion starting from one loop. The MF approximation 
for ChPT implies then the tree level approximation for 
the non-linear sigma model with 5Ch neglected. Such 
an approximation, however, neglects the Haar measure 
from the start. It is therefore hard to expect that such 
an approximation describes correctly the high tempera- 
ture regime where the chiral invariance is supposed to be 
restored. 

The self-consistency of the MF approximation of the 
non-linear sigma model survives with the one parameter- 
ization only. 

The invariance under the chiral transformations ad- 
mits, within the MF approximation, only the parameter- 
ization given by Eg. \IX.Tl[ l. 

The parameterization based on the dilatation of 
gives 6Ch — 0, does not involve the higher orders ChPT 
loops, and allows to work in the continuum limit with 
finite quantities only. 

The effective interaction terms in the Lagrangian 
which appear due to the presence of the Haar measure 
have been discussed earlier in QCD (see fi^l and refer- 
ences therein). 

The SU{2) group has a finite group volume. The in- 
tegration range of u fields is therefore be restricted. Ac- 
cording to the current paradigm, one can extend the in- 
tegrals over LO from — oo to -l-oo within a perturbation 
theory framework. The modification of the result is con- 
nected to the integration over large field fluctuations, so 
the variance has, apparently, a non-perturbative nature 
and does thereby not affect the perturbation series. As a 
matter of fact, this justifies the standard loop expansion 
in ChPT. 



X. SUMMARY 

In this work, we discussed analogy between the second- 
class constraints systems and gauge theories with the 
equivalent structure of gauge generators and gauge-fixing 
conditions. Given the symplectic basis for the constraint 
functions exists globally, the second-class constraints sys- 
tems can be interpreted as gauge invariant systems in the 
unconstrained phase space. Such systems can be quan- 
tized using the methods specific for gauge theories. 

The second-class constraints Qa split in the symplectic 
basis into canonical pairs {xa,^b) satisfying {xa, Xb} ~ 
0, {Q,A,i^B} ~ 0, and {xa,^b} ~ Sab- The constraint 
functions {xa,^b) can be transformed further, as dis- 
cussed in Sects. 4 and 5, to fulfill the Poisson bracket 
relations in the strong sense in an entire neighborhood of 
any give point of the constraint submanifold. The Hamil- 
tonian function can also be modified to be identically in 
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involution with the constraints. The new constraints de- 
fine the same constraint submanifold, whereas the new 
Hamiltonian and its first derivatives coincide with the 
original ones on the constraint submanifold. The con- 
strained dynamics is thus not modified. The new con- 
straint functions XA and ^Ia are interpreted as gauge- 
fixing conditions and first-class constraints associated to 
gauge transformations. 

We do not provide any criterion for what part of the 
constraints Qa describes the gauge-fixing conditions and 
what part describes the first-class constraints associated 
to gauge transformations. By contrary, we argue that 
transition amplitudes of the quantum theory do not de- 
pend on the interpretation of Qa- 

The Dirac's supplementary conditions ilyi^' = de- 
pend on the way the constrains Ga were split. These con- 
ditions are equivalent to Xa"^' = 0, since the correspond- 
ing Wigner functions coincide on the constraint submani- 
fold. The supplementary conditions for the Wigner func- 
tions, furthermore, can be made to be explicitly invariant 
with respect to possible transformations of the constraint 
functions. The ambiguity reflects the freedom in extrap- 
olation of the Wigner functions from the constraint sub- 
manifold into the unconstrained phase space. 

We showed, finally, that the proposed quantization 
scheme applies to an n — 1-dimensional spherical pendu- 
lum, which represents a mechanical version of the 0{n) 
non-linear sigma model. For this model, we demon- 
strated the existence of an underlying gauge symmetry 
which is the dilatation of the coordinates (/>". The con- 
straints appearing within the Hamiltonian framework are 
of the second class. If one starts, however, from an equiv- 
alent Lagrangian in which the underlying gauge symme- 
try is set up explicitly, the same constraints appear as a 
gauge-fixing condition and a constraint associated to the 
dilatation symmetry. It shows that the interpretation of 
second-class constraints is a matter of convention. 

For holonomic systems, the quantization method based 
on construction of the gauged model does requires neither 
auxiliary canonical variables nor extended configuration 
space. 

For second-class constraints systems, the underly- 
ing gauge symmetries induce in the configuration space 
transformations depending on velocities and involve aux- 
iliary variables. The holonomic systems admit the nat- 
ural gauged counterparts in the original configuration 
space. 

The initial configuration space is, in general, too nar- 
row to reflect the gauge invariance of a system de- 
scribed in terms of the generalized Hamiltonian dynam- 
ics. Gauge invariant quantities involving auxiliary vari- 
ables, furthermore, do not belong to the set of physi- 
cal observables, as distinct from the usual gauge theo- 
ries. The equivalence of the first-class constraints sys- 
tems with the ordinary gauge systems is therefore physi- 
cally not complete and restricted to systems of point par- 
ticles under holonomic constraints, non-holonomic con- 
straints satisfying the Frobenius' condition, and systems 



with one primary constraint only. 

The path integral representation for the evolution op- 
erator of the 0{n) non- linear sigma model was con- 
structed in Sect. 9 by solving the constraint equations. 
The equivalence with the quantization methods based 
on the reduction to the equivalent gauge systems was 
demonstrated. 

After finishing this work we got to know about works 
^3iEHISE3 where projection formalism is discussed. 
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APPENDIX A: TWO-CONSTRAINTS FORM OF 
GAUGED SPHERICAL PENDULUM 

Let us consider another example. The system dis- 
cussed in Sect. 3 is equivalent to a system described 
by a vector 0* — ((/>", on the "light-cone" submani- 
fold = (j)0(/)0 _ ^a^a = 0. The constraint = 1 is 
equivalent to the constraint (jP = 1, while the transfor- 
mations ljII.12|) are equivalent to Lorentz boosts of the 
null vector 0* along the vector 



(j)' ^ (1)'' = K'm^P = exp(6')(/)* 



(A.l) 



where the "boost velocity" v = tanh(0). The gauge 
invariance of the £» with respect to the dilatation re- 
flects gauge invariance of the system with respect to the 
Lorentz boosts. 

The "light-cone" Lagrangian £2 can be constructed by 
considering the requirement of the conditional maximum 
of the action 

max{ / £*di}|0o^==i = max{ / £2^0102=0,00=1- 

0",0" J 0^,0* J 

(A.2) 

If the is treated as a gauge parameter, the problem 

simplifies. It is sufficient to require 

max{ / £*cJt} = max{ / £2dt}|02^o- (A.3) 

0",0" J 0^,0^ J 

The £2 can be chosen as a straightforward extension 
of £*: 

£2 = -iG.,(^)0VV(W)' - (A.4) 
where 77 = (1, 0, 0), 



My 



(A.5) 
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and gij = diag(l, -1, -1), 
The tensor Gij obeys 



g^j■ 



ri'G,,{cp) = 0''G,,(0) = 0, (A.6) 

G.,(0)Gi(</>) = GM- (A.7) 

It is invariant with respect to the dilatation (|A.ip and 
the shifts 



+ erf 



(A., 



where e is an arbitrary parameter. 

Lagrangian l|A.4|l is weU defined for cjp' ^ 0. It is in- 
variant with respect to the dilatation (jA.4p and the shifts 



The 0^ and the rjcf) are thus gauge functions. They 
are not fixed by equations of motion and can be selected 
to satisfy admissible constraints. We consider therefore 
Lagrangian ljA.4|l without imposing any constraints. The 
initial and final conditions 0^ = and rjcj) = 1 are gauge- 
fixing conditions. 

The canonical momenta corresponding to the are 
defined by 



8f 

They satisfy the primary constraints 

TTi - GijTT^ W 

which are equivalent to two ones: 
= 07r « 0, 

0,2 — rjTT 0. 



(A.9) 



(A.IO) 



(A.ll) 
(A.12) 



The primary Hamiltonian can be obtained with the use 
of Legendre transformation: 

n = -^((#)' - ^')G''W7T^^T,. (A.13) 

The Poisson bracket for canonical coordinates and mo- 
menta have the form 



(A. 14) 



The primary constraints are stable with respect to the 
time evolution: 



{n,,n} = 0, 
{n2,n} - 0. 



(A.15) 
(A.16) 



The Hamiltonian Ti. is gauge invariant. The primary con- 
straints are of the first class: 



(A.17) 



The generators of gauge transformations constitute an 
algebra. 



The relations 

{<l>\n,} = 4,\ {4,\02}^v\ (A.18) 

{^,,0i] = ~7r„ {^„1]2} = (A.19) 

show that the ili generates the dilatation of the 0* and 
the TTi, while the generates time-like shifts of cj)^ . 
The gauge-fixing conditions 

1 



(A.20) 

X2 = r;0-l (A.21) 
generate the following transformations: 



{0\X2}-O, (A.22) 



K,xi}- (^'^)2t'^t K,X2} = -^^ (A.23) 

They are identically in involution with the Hamiltonian: 



= 0, 

{X2,H} - 0. 



(A.24) 
(A.25) 



The equations of motion generated by the primary 
Hamiltonian look like 

= W,n} = -{{ii(t>? - 4>^)G''{4>>,, (A.26) 

= « -(<^,; - 77,(77(/.))GJ'=(</))7r,7rfe. (A.27) 

The main inference is that starting from the different 
Lagrangian, an equivalent first-class constraints system 
was constructed. There is no principal distinction be- 
tween the system described here and the one discussed 
in Sect. 3. In particular, one can solve the constraints 
X2 = and 1^2 = to remove the canonically conjugate 
pair (0°,7ro) from the Hamiltonian. The system of Sect. 
3 would be reproduced then explicitly. 



APPENDIX B: SYMPLECTIC BASIS FOR 
FIRST-CLASS CONSTRAINTS 

In Sect. 4, an equivalent system of second-class con- 
straints satisfying involution relations Eas. (jIV.23|) in the 
strong sense in an entire neighborhood of a given point 
of the constraint submanifold has been constructed. The 
existence of an equivalent Hamiltonian identically in in- 
volution with the new constraints has also been demon- 
strated Eq. (|IV.24I) . The equivalence means that the con- 
straint submanifolds and the phase space fiows on the 
constraint submanifolds of the dynamical systems coin- 
cide. 

Similar statements for first-class constraints systems 
are prov ed in Refs. 0, [s^l- The arguments of 
Refs.|lOl lU cannot be extended to second-class con- 
straints without additional assumptions. Let us discuss 
the restrictions. 



23 



To replace the constraint functions Qa = 1, ■■■,rn) 
by equivalent constraint functions fi^ which are identi- 
cally in involution {0,a,^b} = 0, one can solve equa- 
tions Qa = with respect to first m canonical momenta 
PA = Pa where Pa are functions of the n canonical coor- 
dinates and the remaining n — rn canonical momenta. To 
have the constraints Qa = resolved, one has to require 



det 



^0. 



(B.l) 



The phase space has a dimension 2n, n > m. The func- 
tions riA — PA — Pa vanish on the submanifold CIa = 
only. The Poisson bracket {^a, ^b} vanishes weakly and 
does not depend on the first m canonical momenta, there- 
fore it vanishes identically. 

The same prescription can be used to construct the 
functions XA- The Poisson bracket {xajXb} van- 
ishes identically. The new constraint functions satisfy 
det{xA,^B} ^ 0. 

One can define an equivalent Hamiltonian 7i. Let us 
substitute pa = Pa to the original Hamiltonian H' of the 



second-class constraints system. The resulting Hamil- 
tonian Ti is first class with respect to the constraints 
Ga = 0, so {flA,ii} « 0. The difference Ti — H' vanishes 
on the constraint submanifold ^a = 0. The Hamiltonian 
H does not depend on the first m canonical momenta, so 
the Poisson bracket does not depend on these canonical 
momenta either. The Ti. is therefore identically in involu- 
tion with the Qa- The similar procedure can be applied 
for the xa- 

We should get finally constraint functions xa and Qa 
identically in involution with the H. 

The above arguments do not apply if some constraint 
functions do not depend on pA {qa)- The bracket 
{(Iat^b} can, e.g., be proportional to xc and xc can 
in turn be independent on pA- In such a case, the weak 
equation {Qa,^b} ~ does not convert into the strong 
one, although the both sides do not depend on pA- The 
similar restrictions appear in the construction of xa and 
7i. The systems under holonomic constraints have, in 
particular, constraints xa = which do not depend on 
the canonical momenta. 
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